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ABSTRACT 

We measure the redshift-space power spectrum P(k) for the recently completed IRAS 
Point Source Catalogue (PSC) redshift survey, which contains 14500 galaxies over 84% 
of the sky with 60 micron flux > 0.6 Jansky. Comparison with simulations shows that 
our estimated errors on P{k) are realistic, and that systematic errors due to the finite 
survey volume are small for wavenumbers k ^ 0.03 /iMpc -1 . At large scales our power 
spectrum is intermediate between those of the earlier QDOT and 1.2 Jansky surveys, 
but with considerably smaller error bars; it falls slightly more steeply to smaller 
scales. We have fitted families of CDM-like models using the Peacock-Dodds formula 
for non-linear evolution; the results are somewhat sensitive to the assumed small-scale 
velocity dispersion cry. Assuming a realistic cry sa 300 km s -1 yields a shape parameter 
r ~ 0.25 and normalisation bas ~ 0.75; if cry is as high as 600 km s -1 then T = 0.5 
is only marginally excluded. There is little evidence for any 'preferred scale' in the 
power spectrum or non-Gaussian behaviour in the distribution of large-scale power. 
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1 INTRODUCTION 

As is well known (e.g. Peebles 1980), the power spectrum of 
the galaxy distribution on large scales is of great importance 
for testing cosmological models, since it can be related to the 
initial conditions by linear perturbation theory. The power 



spectrum is P(k) oc k n with n ~ 1 as suggested by infla- 
tion, then for consistency with the COBE DMR results the 
present-day power spectrum must show a turnover to this 
slope at k ^ 0.02 /iMpc - , close to the largest scales acces- 
sible to current galaxy surveys. There is marginal evidence 
for suc h a turnover in the APM data ( Baugh_fcj]fstjrUii^ 



shift surveys7ToTal3r^Th*e~CrA"Tecrsh7it*^^ 



1994), the QDOT survey (Feldman, Kaiser fc Peacoc k 1993 
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spectru m has been estimated from a variety of galaxy re d- 1994a|; pVladdox, Efstathiou fc Sutherland 1996| ; [Tadros, Ef- 



stathiou fc Dalton 1998) 



hereafter FKP), the Las Campanas red shift survey (Lin et 
19961), and the 1.2 Jy IRAS survey ( [Fisher et al. 1993| ). 



Also, the real-space power spectrum has been inferred from 



the APM Galaxy Su rvey (Baugh fc Efstathiou 1993 Baugh 



& Efstathiou 1994a ) by inversion of both the angular corre- 



lation function and 2-D power spectrum. 

Despite these substantial surveys, there are still consid- 
erable uncertainties in the shape of the power spectrum on 
large scales, since most of these surveys contain only a small 
number of independent structures, while the largest one (Las 
Campanas) has a slice-like geometry which complicates the 
estimation of the power spectrum. If the primordial power 



Also, it is valuable to measure the power spectrum from 
surveys with different selection criteria (e.g. optical & IRAS 
selection) . This is of considerable interest since the observed 
power spectrum is measured from the density field of galax- 
ies, whereas theory predicts the power spectrum of the mass 
distribution. The process of galaxy formation is poorly un- 
derstood, so the observed P g (k) may differ from P m (fc), pos- 
sibly in a complex way; indeed, since it appears that IRAS 
galaxies and optical galaxies have different small-scale cor- 
relation amplitudes, at least one of these cannot trace the 
mass. A simple 'linear bias' model is often assumed, in which 
S g — b8m for some constant 'bias factor' b which may depend 
on galaxy type; this model predicts that P(k) for optical 
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and IRAS galaxies should differ by a multiplicative factor 
of (bo/bi) 2 . Such a model is reasonable since it has been 
shown by several authors (e.g. Fry & Gaztanaga 1993; Cole 
et al. 1998; Mann et al. 1998) that if the galaxy density is a 
(possibly complex and stochastic) function only of the local 
mass density on scales ,$ l/i _1 Mpc, then the effective bias 
parameter defined by b(r) = -^/£ 9 (r)/£ m (r) tends to a con- 
stant on large scales, so such a 'local' bias cannot alter the 
large-scale shape of the power spectrum. 

Another motivation for measuring the power spectrum 
from a large sample of IRAS galaxies is that there appears to 
be a marginal discrepancy between the power spectra from 
the previous QDOT (Fcldm an, Kaiser fc Peac ock 1993, here- 
after FKP) and 1.2 Jansky ( ]Fisher et al. 1993| ) IRAS surveys, 
with the amplitude of P(k) from QDOT being roughly a fac- 
tor of 2 higher at large scales. A counts-in-cells comparison 
of th e two surveys doe s not reveal any obvious systematic er- 
rors (Efstathiou 1995), but it is interesting to check whether 
these differences are consistent with sampling fluctuations in 
one or both surveys. 

In this paper, we estimate the redshift-space powe r 
spectrum from a new redshift survey (Saunders et al. 199C) 
of some 14,500 galaxies over 84% of the sky, selected at 60/im 
from the IRAS Point Source Catalog (PSC). A variety of 
analyses from this survey will appear shortly; the topology 
of the density field has been analysed by Canavezes et al. 
(1998), the correlation function is analysed by Maddox et al. 
(in preparation) , the dipole is estimated by Rowan- Robinson 
et al. (in preparation), a reconstruction of the peculiar veloc- 
ity fiel d is given by Branchini et al. (1998), and the redshift- 



space d stortions are estimated using spherical harmonics by 
Tadros et al. (in preparation). 

The plan of this paper is as follows: in §2 we sum- 
marise the construction and properties of the survey; in §3 
we present the power spectrum estimates, and we compare 
these with results of N-body simulations in §4. In §5 we com- 
pare our results with other surveys and some parametrised 
cosmological models, and also set limits on non-Gaussian 
behaviour and periodicities. 



2 THE PSC REDSHIFT SURVEY 

The construction of the PSC redshift survey (hereafter 
PSCz) is described in detail elsewhere (Saunders et al. 1996; 
Saunders et al. in preparation), but we summarise the main 
points here. The aim of the survey is to obtain redshifts for 
all galaxies with 60/im flux /so > 0.6 Jy over as much of the 
sky as feasible. The starti ng point for the survey is the Q MW 
IRAS Galaxy Catalogue ( Rowan- Robinson et al. 1991), but 



with modifications to extend the sky coverage and improve 
completeness. We have relaxed the IRAS colour criteria for 
galaxy selection, and we have added in additional sources 
in the '2-HCON' sky as follows: the IRAS satellite covered 



most of the sky with 3 hours-confirmed scans (HCONs) ( Be 



ichman et al. 1988), while about 20% of the sky had only 



2 HCONs. Since a source must be detected in 2 separate 
HCONS for inclusion, the PSC catalogue may be less com- 
plete in the 2-HCON regions. Thus, in the 2-HCON sky we 
added sources to our target list which had a 1-HCON detec- 
tion in the 'Point Source Reject' file and also had a matching 
entry in the IRAS Faint Source Catalog. 



These relaxed selection criteria allowed more contam- 
ination of the target list by non-galaxy sources, but these 
were excluded using APM or COSMOS scans of the POSS 
and UKST sky survey plates. If the APM or COSMOS data 
showed no 'obvious' galaxy candidate near the IRAS source, 
we visually inspected the plate and attempted to classify the 
source, rejecting it if it showed an obvious Galactic counter- 
part e.g. an HII region, planetary nebula, dark cloud etc. We 
also exclude very faint galaxies (Bj > 19.5) from the redshift 
survey since measuring their redshifts is time-consuming, 
and they are usually at z > 0.1 and hence have little effect 
on most of the desired analyses. 

The sky coverage of the survey is the whole sky, exclud- 
ing areas with less than 2 HCONs in the IRAS data, regions 
with optical extinction Av > 1.42 mag as estimated from 
the IRAS 100/im maps, and two small areas near the LMC 
and SMC. The resulting coverage is 84% of the whole sky. 
(An extension to 93% sky coverage is in progress, using a 
combination of K-band snapshots and HI redshifts). 

Our 2-D source catalogue contains 17060 IRAS sources 
in the unmasked sky. Of these, 1593 are rejected as objects 
in our own Galaxy (e.g. cirrus, bright stars, reflection nebu- 
lae, planetary nebulae etc), or as multiple entries from very 
nearby galaxies 'broken up' by the IRAS point source de- 
tection scheme. Another 648 sources are rejected either as 
very faint galaxies (~ 400) or as sources without an optical 
identification. This leaves 14819 galaxies in the 'target' list, 
and redshifts are now known for 14539 of these (98%). 

Of these redshifts, ~ 6500 are from a combination of 



the 1.2 Jy s urvey (Fisher et al. 1995) and QDOT (Lawrence 
et al. 1998), and ~ 3000 are from other publications and 



private communications. A further 4115 redshifts were mea- 
sured by us for this survey, using 49 nights at the Isaac New- 
ton Telescope, 18 nights at the Cerro Tololo 1.5-meter, and 
6 nights at the Anglo-Australian Telescope, between 1992 
January and 1995 July. Details of the observations and data 
reduction will be given elsewhere (Saunders et al., in prepa- 
ration). The error on our redshifts is typically 150kms _1 ; 
for the literature redshifts it is somewhat smaller. The me- 
dian redshift of the sample is « 8500 km s -1 , though there 
is a long 'tail' extending to > 30000 kms -1 due to the broad 
luminosity function of IRAS galaxies. 



3 POWER SPECTRUM ESTIMATION 

For the estimates here, we restrict the analysis to the un- 
masked sky with the additional constraint |6| > 10°, since 
the survey may be slightly incomplete below this latitude; 
this gives a coverage of 78% of the full sky. We also set an 
upper redshift limit of cz < 45000 km s" 1 , since the survey is 
incomplete at high redshift as noted above. This gives 13346 
galaxies in the 'default' sample used for the power spectrum 
estimate. 

Since the geometry of the survey is well approximated 
by a sphere, apart from the missing slice near the galac- 
tic plane, we follow the analysis of Feldman, Kaiser & Pea- 
cock (1994, hereafter FKP) with minor modifications. This 
method provides an optimal weighting scheme with redshift 
for estimating the power spectrum of an all-sky survey. More 
sophisticated methods have been suggested by e.g. Tegmark 
(1995); these are very useful for surveys with highly non- 
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spherical geometries but are more complex than necessary 
for our survey. 

We convert the galaxy positions to comoving coordi- 
nates assuming Oo = 1 and redshifts in the Local Group 
frame, and bin the galaxies in a cube of size 950 fc -1 Mpc with 
128 3 cells. The FKP method assigns a redshift-dependent 
weight to each object, 



w(r) = 



>A{\ + P c n(r)) 



(1) 



where n(r) is the mean galaxy density at distance r, P c is 
the estimated power spectrum (at some scale to be deter- 
mined), and A is a normalisation constant (see later). We 
use a parametric fit for the selection function determined us- 
ing the method of Springel & White (1998); here this takes 
the form 



i[z) = n,y (1 - a) /(l + y 1 ) (P/l) 

y = z/z*, z„ = 0.0318 

a = 1.769, 13 = 4.531, 

n, = 8.76 x 10" 3 h 3 Mpc~ 3 ; 



(2) 



7 = 1.335 



these values are appropriate for / > 0.6 Jy; for other flux 
limits we simply scale z* by */(0.6//n m )- 

We have made three refinements to the FKP estimator: 



firstly, we define the ratio of densities of real and random 
catalogues a' = ^2 g w i/^2 s w j, where Wi is the weight of 
the ith object and t he sums run over galaxies and random 
points respectively (Tadros & Efstathiou 1996), instead of 
a — Ng/N r as in FKP (where N g ,N r are the numbers of 
galaxies and randoms respectively). Secondly, we compute 
the shot noise using 



-Pshot = ^ Wi + a' 2 Wj, 



(3) 



where the two terms are the contributions from galaxies and 
random points respectively. The shot noise definition in FKP 
Eq. 2.4.5 was P s hot = a(l + a) ^ wi] there the first-order 
term in a is the 'expected' shot noise from the galaxies given 
many realisations of the given selection function, while the 
first term in our definition is the 'actual' shot noise in the 
data. This makes negligible difference at large scales, but we 
find from simulations that Eq. ^ gives substantially smaller 
errors in the estimated power spectrum at small scales (large 
fc), because the shot noise term is substantial here and the 
'actual' shot noise from the galaxies may differ significantly 
from its expectation value estimated from the selection func- 
tion. 

The third refinement is that we use a normalisation con- 
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Figure 2. Weight and window functions vs. rcdshift for P c = 
8000 h _3 Mpc 3 , in arbitrary units. The solid line shows the weight 
w(z), and the dotted line shows the real-space window function 
n(z)w(z). The dashed line shows n(z)w(z)dV/dz, i.e. the 'effective 
volume' of the survey per unit redshift, and the dot-dash line is 
the integrated effective volume below redshift z. 



imum redshift of 30000 km s~ , or a galactic latitude cut of 
20°, and changing the flux threshold. Most of these changes 
have a negligible effect on the results, with the exception of 
changing the flux limit: the subsample with feo > 0.8 Jy has 
a slightly lower amplitude of P(k) at all scales; the difference 
is mostly within the individual la error bars but is evident 
over a range of k. The slice with 0.6 < fm < 0.8 Jy has 
correspondingly higher amplitude. This might be suggestive 
of systematic errors in the catalogue near the flux limit, as 
suggested by Hamilton (1996); however we have investigated 
the correlation function £(cr, n) as a function of projected 
and redshift separation and find negligible evidence for a 
"Hamilton effect" i.e. elongation of the correlation function 
out to large redshift separation. The correlation functions of 
the subsamples above and below 0.8 Jy show a similar differ- 
ence in amplitudes to the P(k)'s, but the difference is most 
pronounced at small to intermediate scales, rather than at 
large scales as might be expected from systematic errors; 
thus this effect may be a sampling fluctuation or possibly a 
dependence of clustering amplitude on intrinsic galaxy lu- 
minosity. Hereafter we use the full sample to 0.6 Jy with 
the caveat that the reason for this slight flux-dependence is 
not yet clearly understoo d. T he effect on the derived cos- 
mological parameters in 



5.2 



is similar to or less than the 
uncertainties arising from the small-scale velocity dispersion 
etc. 



vention given by Eq. Ac ; see Appendix A for a discussion of 
the normalisation. 

The estimated power spectrum P(k) is then given as in 
FKP, by 



F(r) = w(r)[n g (r) - a'n s (r)] 



F(k) 
P(k) 



<frF(r)e 



(4) 
(5) 
(6) 



where n g ,n s are the number densities of galaxies and ran- 
dom points in cubical cells, and P(k) is just the unweighted 
average of P(k) over a spherical shell with mean radius k. 

The optimal weighting scheme depends on the ac- 
tual value of P(k), so the procedure is slightly cir- 
cular in principle. We have used values of P c = 
2000,4000,8000, 16000 /i~ 3 Mpc 3 ; estimates of the power 
spectrum for each value of P e are shown in Figure |l|. 

We see that changing P c changes the size of the error 
bars, but there is little systematic difference in the resulting 
estimates of P{k). This is as expected since FKP showed 
that any choice of P c gives an unbiased estimate of P(k) 
(apart from the convolution effects at small k discussed be- 
low), but just weights different redshift shells differently - 
larger P e gives relatively more weight to more distant shells. 
We adopt P e = 8000 /i" 3 Mpc 3 as the default value for the 
remainder of the paper. The resulting weights are illustrated 
in Figure ^| the solid line shows the weight function w(z), 
and the dotted line shows the real-space window function 
n(z)w(z). Also shown are the differential and cumulative 
contributions to the survey 'effective volume' per unit red- 
shift. 

We have explored varying many of the selection criteria, 
e.g. varying the model of the selection function, using a max- 



3.1 Observed vs True P(k) 

As is well known, the finite size of the survey volume means 
that the estimated power spectrum is a convolution of the 
true power spectrum with the squared Fourier transform of 
the real-space window function; e.g. Eqs. 2.1.6 and 2.1.10 of 
FKP give 



(P(k)) = (2tt)- 3 / d 3 fc'P(k')|G(k-k') 



G(r) = n(r)w(r) I / d 3 rn 2 (r)w 2 (r) 



and the normalisation is defined so that 



d 3 fc|G(k)| 2 = (2tt) 3 



(7) 



-1/2 



(8) 



by Parseval's theorem. This convolution is a significant prob- 
lem for slice-like survey geometries with a highly anisotropic 
window function; but since our survey is large in all 3 di- 
mensions, the window function is narrow. For our standard 
weighting with P e = 8000 /i~ 3 Mpc 3 , the window function 
|G(k)| 2 is illustrated in Figure H for three axes in Galac- 
tic coordinates, and the angle- average. The window func- 
tion is roughly approximated by a Gaussian exp(— k 2 /2k^) 
with ko ~ 0.006 ftMpc" 1 at small k, with a roughly fc -4 tail 
arising from the survey mask. Therefore, the effect of the 
convolution on our estimates is small except at the largest 
scales k £ 0.03 /iMpc -1 . This is illustrated in Figure [I| the 
solid lines show the fractional contribution to the measured 
P(k) as a function of 'true' wavenumber k' , i.e. Eq. |^ av- 
eraged over directions of k, k' for five values of observed 
k = 0.033,0.066, 0.1, 0.2, 0.3 h Mpc~\ assuming a CDM-like 
model P(k') with the parameter F = 0.3 (see Eq. |ll] be- 
low). We have not attempted a deconvolution here, since 
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k (h Mpc 

Figure 3. The k-space window function |G(k)| 2 . Points show this quantity for k parallel to the Galactic x, y, z directions (as labelled). 
The solid line shows the direction-average over spherical shells of radius k. The dashed line shows a Gaussian exp(— fc 2 /2fcjj) with 
fco = 0.007 h Mpc -1 (this is for illustration and is not a fit). 
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Figure 4. The effect of convolution, i.e. the contribution to 'observed' power P(k) from 'true' wavenumber k' in Eq. Q summed over 
directions of k, k', per unit lnfc'. Values are shown for observed k = 0.033, 0.066, 0.1, 0.2, 0.3 feMpc -1 (dashed lines). The y-scale is 
arbitrary. 



the convolution effect is only important at small k where 
the estimates are becoming noisy. 

Another effect which causes the measured P(k) to devi- 
ate systematically from its true value is that the mean den- 
sity of galaxies is not known independently but is estimated 
from the survey. This leads to the constraint P(0) = 0, and 
the convolution above means that P(k) will also be under- 
estimated for small but non-zero k. This effect was noted by 
Peacock & Nicholson (1991), and has been evaluated ana- 
lytically by Tadros & Efstathiou (1996) for the special case 



of a volume-limited survey; their Eq. A4.2 gives 



<P(k) - P(k)) 



V 



\W(k) 2 



(9) 



5]|w(k')! 2 P(k') 

. k' 



|W(k) 2 



For unequal weights as here, the expression is complex and 
best evaluated numerically; we have computed this for two 
models of P(fe); a T = 0.3 CDM model, and an n = -1.2 
power-law, with a bend to n — at k < 0.01 /iMpc -1 . For 
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Figure 5. The effect of convolution and the normalisation con- 
dition P(k) = on two model power spectra (solid lines); the 
upper solid line shows a power law P(k) <x fc -12 with a break 
at k = 0.01 hMpc" 1 ; the lower solid line shows V = 0.3 CDM. 
The dashed and dot-dash lines show the mean recovered power 
spectra for the two cases. Points show the measured values for 
PSCz. 



each model, we generate a Gaussian random density field 
with the assumed P(k), multiply by the PSCz window func- 
tion n(r)w(r), set the mean of the windowed density field 
to zero, and compute the power spectrum of the resulting 
density field using the given w(r). The mean of the recov- 
ered power spectra from 30 realisations of each model is 
shown in Figure ^[ We see that the bias is serious only for 
k £ 0.02/i Mpc" 1 ; note in particular that for the power-law 
model, the recovered power spectrum lies well above the 
PSCz measurements at k ~ 0.01 - 0.03 h Mpc" 1 . The data 
points at k ~ 0.013and0.02 h Mpc" 1 are noteworthy here; 
although their error bars appear large due to the log scale, 
none of the 30 realisations of the power-law model gave P(k) 
as low as the PSCz data at these k. 

Thus, we can be confident that the flattening of the 
observed power spectrum below k ~ 0.06 h Mpc" 1 is a real 
feature, not an artefact of the finite volume or normalisation. 

We compute the covariance matrix of P(k) using 
Eq 2.5.2 of FKP. In practice, it is infeasible to evaluate this 
directly since it contains ~ N 6 terms where N = 128; thus in 
practice, we assign each of the N 3 wavevectors to its bin in 
fc; for each k, k' we pick ~ 10 6 random pairs of wavevectors 
k, k' in the appropriate bins, and evaluate the sum accord- 
ingly. It is found that ~ 10 6 pairs are necessary to ensure 
that the random errors are small, otherwise the resulting 
matrix can become non-positive-definite. 

A final effect on P(k) is the 'binning factor' noted by 
Baugh & Efstathiou (1994b), which causes an underestimate 
of small-scale power due to the galaxies being binned into 
finite-size cells before the Fourier transform. This 'smooths' 
the observed density field over the bin size. The size of this 
effect depends on the slope of the true power spectrum; a 
second-order approximation is given by Peacock & Dodds 



(1996) Eq. 20, and is P(k)/P(k) w [l + (kl) 2 / '12] \ where 
I is the size of the unit cells in the FFT; this becomes inac- 
curate for kl ^ 2. We find that for n~-la better approx- 
imation useful for kl < n is 

P(k) /P(k) w [l + (kl f /12 - (fc/) 4 /220] _1 . (10) 

Since this correction is slightly model-dependent, we 
have not applied it to the estimates in Figure [jjor Figure [| 
but we have applied it before the model fitting in the fol- 
lowing Section. 



4 COMPARISON WITH SIMULATIONS 

As a check of the code, and to assess whether the resulting 
error bars are realistic, we have generated simulated 'PSCz' 
surveys from large N-body simulations from 3 cosmological 
models and computed their power spectra as above. The 3 
models are 'standard' CDM (SCDM), CDM with a cosmo- 
logical constant (ACDM), and a mixed cold + hot dar k mat- 



ter model (MDM): the simulations use a P3M code (Croft 



fc Efstathiou 1994 ) and their parameters are listed in Ta- 
ble y. We place an 'observer' at a random location in the 
cube, wrap the simulation using periodic boundary condi- 
tions, and then select 'galaxies' in redshift space as random 
particles using the selection function of Eq. |2j and a 'mask' 
of the same shape as the real PSCz mask. We generated 
a total of 27 simulated 'PSCz' surveys for each model; for 
each model we used 9 different runs of the P3M code, and 3 
different observer locations for each run. 

Figure ^| shows the power spectra of the full box for 
each simulation, compared with results from 9 simulated 
PSCz surveys. As expected, our estimator P(k) recov- 
ers the 'true' redshift-space power spectrum quite well for 
k ^ 0.02 h Mpc" 1 ; the simulated surveys slightly underes- 
timate the true redshift-space P(k) on intermediate scales 
due to the convolution with the survey window function, but 
the effect is small. Figure |^ shows the mean of the FKP er- 
ror bars from 9 simulations (triangles), compared with the 
'real' uncertainty estimated from the rms scatter between 
the 9 simulations (crosses). Clearly the FKP error estimates 
are a reasonable approximation to the 'real' errors in P(k), 
though they appear to underestimate the actual errors by 
~ 20%. 

In Figure ^ we compare the data to the mean of 27 sim- 
ulated PSCz surveys for each model. We find that all three 
models give a reasonable match to the shape of the observed 
power spectrum; the ACDM model has somewhat too high 
an amplitude and would require an antibias b < 1. However, 
this could be remedied by lowering h somewhat and raising 
flo, keeping T constant; this would reduce the implied ag for 
COBE normalisation, and improve the fit. Alternatively, a 
modest degree of 'tilt' with primordial spectral index n < 1 
would similarly reduce the COBE-normalised ag. 

Although the SCDM model has much less large-scale 
power than the others in real space, the high COBE normal- 
isation gives two effects: an enhancement of power on large 
scales by a factor « 1.86 from the Kaiser (1987) redshift- 
space distortion, and a suppression of small-scale power 
from the resulting large peculiar velocities. These two ef- 
fects combine to bring the redshift-space P(k) of this model 
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Table 1. Estimated power spectrum 



k 

(ftMpc" 1 ) 


P(k) 
(h" 3 Mpc 3 ) 


Error 


k 

(hMpc" 1 ) 


P(k) 
(/i" 3 Mpc 3 ) 


Error 


0.0066 


16110 


21200 


0.0839 


8690 


1165 


0.0132 


3981 


5774 


0.1056 


4548 


625 


0.0198 


3503 


4039 


0.1329 


3786 


425 


0.0265 


13514 


5960 


0.1674 


2591 


285 


0.0334 


14189 


4861 


0.2107 


1466 


193 


0.0420 


12459 


3428 


0.2653 


1122 


151 


0.0529 


10241 


2321 


0.3339 


748 


119 


0.0666 


8059 


1488 


0.4204 


545 


89 



This table shows the estimated power spectrum P(k) with weight P c = 8000 /i~ 3 Mpc 3 , and associated la errors. These data points 
have been 'corrected' for finite-size bins using Eq. |l(j with I = 950/128 h~ 1 Mpc. Covariances between consecutive entries are substantial 

for small k but negligible for k ^ 0.1/iMpc -1 . 



Table 2. N-body Simulation Parameters 



Name 


N 


Box (r'Mpc) 


^CDM 


f^HDM 


n A 


h 




SCDM 


160 3 


600 


1 








0.5 


1.0 


ACDM 


160 3 


600 


0.2 





0.8 


1.0 


1.0 


MDM 


100 3 


300 


0.7 


0.3 





0.5 


0.67 



into rather good agreement with the data. These effects 
for COBE-normalised SCDM have been previously noted 
by various authors, e.g. Bahcall et al. (1993), although of 
course t his model has serious pr oblems with cluster abun- 
dances (Eke, Cole fc Frenk 1996), large-separation gravita- 



tional lenses (Cen et al. 1994) etc 



5 DISCUSSION 

5.1 Comparison with Other Surveys 

Our observed redshift-space power spectrum for P c = 
8000 /i _3 Mpc 3 is compared with a number of previous mea- 
surements in Figure ^. These come from various catalogues, 
both optical and IRAS selected; all are redshift-space power 
spectra except for the APM data which come from an in- 
version of the 2-D power spectrum. We see that on large 
scales our measurements are well within the range of previ- 
ous surveys, but on intermediate scales ~ 0.2/iMpc" 1 our 
measurements are slightly steeper than the others, notably 
the combined QDOT and 1.2 Jy surveys. This is the cause 
of the low values of T ~ 0.2 for best-fitting CDM-like mod- 
els seen in the next section. It appears that optical galaxies 
have a somewhat higher power spectrum amplitude on inter- 
mediate scales, as expected from the fact that IR-selected 
surveys contain mainly late-type galaxies, and from their 
smaller correlation length ro , but it remains unclear whether 
this persists to large scales; a direct comparison with APM is 
complicated by the redshift-space distortion, and the inter- 
pretation of LCRS is somewhat complicated by the inversion 
from the 2-D to 3-D power spectrum due to the slice-like ge- 
ometry. Future large optical surveys such as 2dF and Sloan 
should greatly clarify this question. 



5.2 Fits to the power spectrum 

In addition to the direct comparison with simulations, it 
is interesting to extract best-fitting values for parametrised 
models of the power spectrum; we use firstly linear theory 
for simplicity, and later the fitting formulae of Peacock & 
Dodds (1996) which account both for non- linear evolution 
of clustering, and the effect of distortions between real space 
and redshift space. 

We use CDM-like models with the initial power spec- 
trum parametrised by F as in Eq. 7 of Efstathiou, Bond & 
White (1993, hereafter EBW), which is 



P(k) 



Bk 



(1 + [afc + («c) 3/2 + {ck) 2 Yf /v 
6.4/rft _1 Mpc, b = 3.0/r/i _1 Mpc, 
1.7/r/j _1 Mpc i/=1.13 



(11) 



In linear theory there are only 2 free parameters: the shape 
parameter T, and the normalisation which may be taken as 
bos, where b is the bias parameter and erg is the rms mass 
fluctuation in an 8/i _1 Mpc top-hat sphere. The % 2 contours 
using the full covariance matrix are shown in Figure O; the 
best fit values are F = 0.19 ± 0.03 and ba 8 = 0.80 ± 0.02. 
These compare to T = 0.19 ± 0.06, ba s = 0.87 ± 0.07 for the 
QDOT survey, as given by FKP Eq. 4.3.3. 

For the non-linear Peacock-Dodds formula, we need a 
total of 6 parameters to specify the present-day redshift 
space power spectrum: the initial mass power spectrum is 
specified as above by F and erg (where as is defined as 
the initial rms mass fluctuation, multiplied by the linear- 
theory growth factor to the present day). The subsequent 
non-linear evolution depends also on f2o and (weakly) on 
Qa- The transformation from real to redshift space depends 
on flo and the bias b, and also on the pairwise peculiar ve- 
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Figure 7. Error estimates for simulated PSCz surveys, with the observed mask and selection function, for the three models in Tabic Q 
Triangles show the mean of the FKP error estimate 6P(k) from 9 simulated surveys; crosses show the 'true' error in P(k) estimated from 
the scatter in the 9 simulated surveys. 




locity dispersion; in the PD formula this is approximated by 
assuming galaxy ve locities a re independent Gaussians with 
dispersion ay = y v\ 2 (r) /2 for some suitable scale r. We 
assume simple 'linear bias' so the galaxy power spectrum 
is b 2 x the matter power spectrum. There is clearly insuffi- 
cient information in the P(k) data to fit all 6 parameters 
separately, so we restrict the parameter space as follows; 



(a) To set fi, we consider either Einstein-de Sitter mod- 
els, with Q = 1, £Ia = 0, treating T as a free param- 
eter (which is a reasonable approximation to e.g. mixed 
dark matter models); or we consider ACDM models with 
Q — T/0.66, Qa = 1 — fi (i.e. assuming F = flh with a 
Hubble constant h — 0.66, consistent with most recent mea- 
surements) ; 
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Figure 9. Comparison of PSCz P(k) with other measured power spectra. The data points as labelled on the figure are: IRAS 1.2 
Jansky from Fisher et al. (1993), QDOT from FKP, QDOT + 1.2 Jansky from Tadros & Efstathiou (1995), Stromlo-APM from Tadros 
& Efstathiou (1996), APM (real-space, deconvolved) from Baugh & Efstathiou (1994), and Las Campanas (deconvolved) from Lin et 
ai.(1996). For clarity, error bars are only shown on PSCz and APM (others are larger), and Stromlo, QDOT and QDOT+1.2 Jansky 
data have been rebinned. Dotted lines show linear-theory CDM with T = 0.2 and 0.5 with erg = 0.8. 



Table 3. Fits of CDM-like models 



n 


Norm (erg) 


cry ( km s *) 


r 


b 


<?8 


Plot? 


Linear 






0.19 





80 


* 


1 


COBE (0.45) 


300 


0.20 





67 




1 


COBE (0.49) 


MJB (368) 


0.22 





68 




1 


COBE (1.00) 


600 


0.40 





71 




1 


Clus (0.52) 


300 


0.20 





66 




1 


Clus (0.52) 


MJB (390) 


0.23 





69 




1 


Clus (0.52) 


600 


0.34 





81 




T/0.66 


COBE (0.97) 


300 


0.16 





66 




T/0.66 


COBE (0.97) 


MJB (457) 


0.16 





70 




T/0.66 


COBE (1.38) 


600 


0.34 





73 




r/0.66 


Clus (0.91) 


300 


0.20 





67 


* 


r/0.66 


Clus (0.81) 


MJB (438) 


0.25 





74 


* 


r/0.66 


Clus (0.73) 


600 


0.31 





83 


* 



Fits to the observed P(k) in the (T, bag) plane. The first line shows linear theory, the remainder use the non-linear PD formula, using 
various choices for setting the other input parameters (f2, Qjy, erg, ay) as a function of T. Column 1: either Q = 1,S1a = or 
Q = 0.66/r, Qj^ = 1 — f2; Column 2: erg defined either by COBE or cluster normalisation; Column 3: cry either fixed to 300 or 
600 kms" 1 , or defined as a function of the other parameters using the MJB model. For the erg and ay columns, the values in brackets 
show the derived values at the best-fitting Y . The * in the last column flags the models plotted in Fig s. and Fll 



© 0000 RAS, MNRAS 000, 000-000 



10 W. Sutherland et al. 




- - Real space 
— z space 
o Avg of sims 
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■ ■ Real space 
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Figure 6. Power spectra for the 3 models in Table ^. For each 
model, the dotted and solid lines show the real and redshift space 
power spectrum of the full simulation box (mean of 9 runs, error 
bars negligible). The circles show the mean of the estimated power 
spectra from 9 simulated 'PSCz' surveys, with lc error on the 
mean. 



(b) To set as, we u se either the cluster nor malisation 
a 8 = 0.52n ( -°- 52+013n) ( [Eke, Cole fe Frenk 1996| ), or we use 
COBE normalisation where as is a function of F using Eq. 6 
of EBW and Q rms = 17 fj,K . 



(c) To set av, we fix either ay = 300 or 600 km s , or 
predict av as a function of V, as using the fitting formula in 
Eqs. 40a, 40b of Mo, Jing & Borner (1997, hereafter MJB). 
A value as hi gh as 600 km s -1 is probably disfav oured by 



observations (Landy, Szalay & Broadhurst 199S ) , but we 



include this as a conservative upper limit to show the effect 
on the derived parameters. 

Having made one choice from each of (a),(b),(c) above, 
this defines Q, Qa, as, ay as a function of T; we then treat 
r and b as free parameters, and fit to the observed P(k) 
data. We choose to use only the data points in the range 
0.021 < k < 0.3 /i Mpc -1 in the fits, since points at lower k 
may be affected by the convolution, and those at higher k are 
subject to large and somewhat uncertain corrections both 
for the peculiar velocity term and the binning correction of 
Eq. 0. 

The best-fitting values of T, bag for each of the above 




k (h Mpc ] ) 

Figure 10. Fits to the measured powe r sp ectrum (points) using 
the Peacock-Dodds formula as in Sec. |5.2| , Vertical dotted lines 
denote the range of k used in the fits. Lines show the fitted power 
spectra as labelled, one for linear theory and three for different 
values of ay, for the case of a low-density universe with SI = 
r/0.66, and cluster normalisation. 



parameter choices are shown in Table q. (We present the fit 
results in terms of (F, bas) rather than (r, b) because 6 and 
as are degenerate in the linear regime, and also for simplicity 
since the contours of equal \ 2 are roughly parallel to the axes 
in the (F, bas) plane.) For the low-Si cluster-normalised case, 
contours of goodness of fit are shown in Figure [ll] for each 
of the three ay assumptions. The derived power spectra for 
the best fits in the same cases are shown in Figure [j^ For 
most of the fits, a fairly small value of F ~ 0.2 is favoured, 
and F ~ 0.5 is quite strongly ruled out; for the cases with 
ay = 600kms _1 , higher values F ~ 0.35 are favoured and 
r = 0.5 is only marginally ruled out, though values of ay as 
high as this are definitely not favoured observationally. 

There is rather little difference between the goodness- 
of-fit for the various choices, though the fits for ay = 
600 km s -1 are somewhat worse. For a given parameter 
choice, the la random errors are typically ±15% in F and 
±4% in bas- In general we see that the systematic uncer- 
tainties from the different choices of fi, erg, ay dominate the 
random errors arising from the error bars on P; the largest 
source of uncertainty is that arising from ay. As expected, 
if we increase ay the predicted small-scale power decreases 
at a fixed F, and thus the best-fitting values of F and bas 
increase to compensate. 

For ACDM models, the value of F at which the COBE 
and cluster normalisations agree is F = 0.17, as ~ 1.0; this 
is interestingly close to our best-fit values of F ~ 0.2. Our fit 
values of bas thus imply b ~ 0.75: this is a significant amount 
of "antibias", but may be plausibly accounted for by the 
deficiency of IRAS galaxies in rich clusters. Thus a ACDM 
model is attractive in that it can simultaneously satisfy three 
constraints (COBE, cluster abundance, P(k)) with one free 
parameter F, if optical galaxies are approximately unbiased 
and IRAS galaxies mildly antibiased relative to the mass. 
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Figure 11. Contours of \ 2 m the (I\ bag) plane, of fi ts to the observed power spectrum. The upper left panel shows linear theory. Other 



panels use the Peacock-Dodds formula as in Sec. 5.2, with f! = T/0.66, cluster normalisation, and different choices of random velocities 



cry as labelled. The cross denotes the minimum \ A , contours are at x 



1,4,9,16,25. 



5.3 Periodicities and Spikes 

There have been a number of suggestions of a 'preferred 
scale' for large-scale clustering, notably by Broadhurst et al. 
(1990), Landy et al. (1996), and Einasto et al. (1997). These 
effects, if real, could arise from a 'spike' in the power spec- 
trum, such as may arise from a baryon isocurvature model 
or non-standard inflation models; or from non-Gaussian ini- 
tial conditions, which could lead to one particular direction 
showing a value of |5(k)| 2 much larger than its expectation 
value P(k); or even from an intrinsic 'preferred direction' in 
the Universe. 

In our data, there is marginal evidence for a 'spike', per- 
haps better described as a 'step', in the power spectrum near 
k ~ 0.08 h Mpc -1 , but this is only about a 2a effect above a 
smooth CDM-like fit, and the scale is significantly different 
from that k ~ 0.06/iMpc" 1 suggested by Broadhurst et al. 
(1990) and Landy et al. (1996); also, we find that changing 
our flux limit to fso > 0.7 Jy causes a substantial drop in 
P(k) at this point, while leaving other points virtually un- 



changed. Thus, we suspect this point may be a statistical 
fluctuation, and there is no conclusive evidence for a feature 
in our power spectrum. From inspection of results from a 
number of realisations of N-body simulations, we find that 
such features quite commonly arise from statistical fluctua- 
tions. 

To test for non-Gaussian effects, we have examined the 
histogram of the ratio of observed to mean power 
|F(k)| 2 /(P(fc) + Pshot) for each wavenumber; for Gaussian 
clustering, this should follow an exponential distribution 
with unit mean. Results for several ranges of wavenumber 
are shown in Figure |l2| and closely follow the exponential 
distribution. The presence of the shot noise somewhat weak- 
ens this test, since for smaller scales P s hot ~ P(k), but at 
large scales this is a sensitive test for non-Gaussian initial 
conditions with a 'tail' to high values. The distributions are 
well fitted by the exponential, so there is no evidence for 
non-Gaussian initial conditions or any 'preferred direction' 
in our survey. 

This also strongly constrains any 'preferred direction', 
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Power / Mean Power / Mean 



Figure 12. Solid lines show histograms of 'observed' power P(k) + P s hot divided by the mean for each k. The dashed line shows an 
exponential distribution with unit mean. The four panels show various ranges in k as labelled. 



as follows. If there exists a strict plane- wave periodicity in 
the universe with dimensionless density contrast a along 
wavevector ko, the true power spectrum contains delta func- 
tions at ±k , P(k) = C[<S(k-k )+<$(k+ko)]. Requiring the 
variance a 2 = a 2 /2 = (2ir)~ 3 J d 3 fcP(k) gives the constant 
C = (2n) 3 a 2 /4. Then, after convolution with the survey win- 
dow function as in Eq. [?], this contributes to the observed 
power spectrum a term of size 



fW(k) = 5L (| G (k - k () )| 2 + |G(k + k )| 2 ) 
Recall that 



\G(k = 0)\ 2 = VJ 



J d 3 rn 2 (r)w 2 (r) 



(12) 



(13) 



for a volume-limited survey with constant n, w this just 
reduces to the survey volume V, so it may be thought 
of as the survey 'effective volume'. For our survey with 
weight function given by P c = 8000 /i~ 3 Mpc 3 , we have 
|G(0)| 2 = 6.2 x 10 7 h~ 3 Mpc 3 : thus even a small amplitude 
periodic wave of a = 0.15 would lead to a large spike in our 
measured P(ko) ~ 3 x 10 5 /i _3 Mpc 3 , well outside the expo- 
nential tail of measured values. We conclude that there is no 



strictly periodic plane-wave structure in our survey volume 
with amplitude larger than 15%. 



6 CONCLUSIONS 

Our conclusions may be summarised as follows: 

(i) The redshift-space power spectrum of the PSCz sur- 
vey is intermediate between those of the earlier QDOT and 
1.2 Jansky IRAS surveys on large scales, though it is slightly 
steeper on small scales. It is stable against variations in the 
galactic cuts and redshift limits etc, though the amplitude 
decreases slightly for a flux cut / £ 0.8 Jy. 

(ii) There is convincing evidence for curvature in the 
power spectrum; the slope changes from the small-scale 
power law n w —1.4 to n ~ on scales k ,$ 0.07 h Mpc -1 . 
This is not an artefact of the finite sample volume or the 
estimation of the mean density from the survey. 

(iii) The best-fitting CDM-like models have F ~ 
0.25, bag ~ 0.7. The uncertainties in these values are mainly 
due to uncertainties in the small-scale velocity dispersion, 
the value of Q etc, rather than statistical errors. 

(iv) There is little evidence for a 'spike' in the power 
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spectrum, and no evidence for large-scale periodicity or non- 
Gaussianity. 
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APPENDIX A: NORMALISATION OF P(K) 

We note an issue concerning the overall normalisation of 
P(k). FKP set a normalisation of their weight function via 
Eq. 2.4.1. It is convenient in practice to set weights via Eq. ^ 
with A — 1, and then later divide all power spectra by a 
constant A. There are several ways of doing this, which give 
similar but not identical results: The LHS of FKP Eq. 2.4.1 
gives 



Ai 



d r n 2 (r)w 2 (r); 



the RHS gives 

A 2 = fty]n(r,)iii 2 (r s ), 



(Al) 



(A2) 



where the sum runs over random points, and a = N g /N r as 
before. Another possible definition is 



A3 = 



(A3) 



where the first sum runs over the cells used in the Fourier 
transfo rm, a is the sum of weights of all random points in 



the ith cell and v is the volume of a unit cell. (This results 
from estimating the window function from the FFT of the 
random points via 



|G(k)| 



A 3 V 



Eik.r, 




(A4) 



requiring ^ k |G(k)| 2 = 1 leads to Eq. |A3 



The relationship between A\,A 2 ,A3 is as follows: the 
number of random points N r is arbitrarily fixed, and we 
define j3 = J d 3 rn(r)/N r to be the ratio of the expected 
number of galaxies from the selection function to the number 
of randoms, which should be similar but not identical to a. 
Thus the expected number of randoms in cell i is just n;t)//3 
where n% = n(ri). For small cells, all randoms in a given cell 
will have equal weight Wi, thus c,: = Wi x Poisson(niv//3). 
Thus 



(A 2 ) = a^^(ci)mw 2 

cells 

= a^^n 2 vw 2 /f3 



(A5) 



Since d is a Poisson variable multiplied by Wi, (c 2 ) = 



Wi [(riiv/P) + nv/ 13], thus Eq. A3 gives 



cells cells / 



P 2 



OTV \ -> 2—2 



(A6) 



Thus we see that the three definitions of A differ by pow- 
ers of a/ P, which is just the ratio of 'observed' to 'expected' 
number of galaxies. Which is more 'correct' is largely a mat- 
ter of choice, but the definition A3 is convenient in practice 
since it leaves the results unchanged if we rescale the selec- 
tion function by a constant n —* en and rescale the weights 
by P c -> PJc. 
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